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NMR measurements of the 17 nuclear quadrupole interactions, 11 vq, show that the charge 
environment of the S=l/2 CU2O3 two-leg spin-ladder layer of hole-doped A14CU24O41 changes dra- 
matically above a certain temperature, T*. This temperature, T*, decreases with additional doping 
and is correlated with the magnetic crossover from the spin gap regime to the paramagnetic regime. 
We demonstrate that these changes in 1 vq are consistent with an increase in the effective hole 
concentration in the CU2O3 two- leg ladder above T*. The effective hole concentration increases 
primarily in the oxygen 2pcr orbitals. 

PACS numbers: 76.60, 74.25.Nf, 74.72 



I. INTRODUCTION 

The interplay between quantum spin fluctuations and 
doped holes in the S=l/2 CuC>2 square-lattice has been a 
major controversy in high T c superconductivity. The dis- 
coveries of the spin-liquid ground state/with spin gap A in 
the S=l/2 CU2O3 two-leg spin ladder p and superconduc- 
tivity under high pressure in hole doped ladders,cl raise 
the same question in a reduced dimension. In this article, 
we report a detailed study of the 17 O and 63 Cu nuclear 
quadrupole interaction tensor, 17,63 z/q, in the CU2O3 two- 
leg ladder layer of the hole-doped ladder-chain compound 
A14CU24O41 (A=La 6 Ca 8 , Sru, SrnCa 3 , Sr 6 Cag). 

The nuclear quadrupole interaction tensor, 17 ' 63 fQ, 
measures the electric field gradient at the position of the 
nucleus, and thus probes the local charge environment. 
We define vq based on the nuclear quadrupole Hamilto- 
nian asu 

n Q = v z Q Wl Hi + 1) + \n{i% + 1 2 )} (ia) 



17 0). eq is the electric field gradient at the site of the 
nucleus along the principle axis z of the tensor, 

^ = E^- ( 2 ) 

where the summation is taken over all charges including 
ions, electrons, etc. The asymmetry parameter, 77, is the 
deviation of the electric field gradient tensor from axial 
symmetry. 

The electric field gradient, eq, is sensitive to changes 
in the amount and symmetry of the charge distribution. 
We demonstrate that the charge environment of the hole- 
doped ladder changes dramatically above a certain tem- 
perature, T*, somewhat below the spin gap (T* < A). 
Our calculations indicate that extra holes moving primar- 
ily into oxygen 2per orbitals are necessary to reproduce 
the temperature dependence of 17,63 ^q above T*. 



II. EXPERIMENT AND RESULTS 
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where Q is the quadrupole moment of the nucleus and I 
is the spin of the nucleus (1=3/2 for 63 - 65 Cu, 1=5/2 for 



We grew single crystal samples of undoped 
La6CasCu2404i (valence of Cu is +2, S=l/2) and 
Sri4Cu24C>4i (amount of holes Pl ~ 0.06 per ladder Cu 
at room temperatureu) using the floating zone technique. 
For SriiCa 3 Cu240 4 i [P L - 0.12) and Sr 6 Ca 8 Cu2 4 4 i 
(Pl ~ 0.17), we uniaxially aligned ceramic powder in 
epoxy along the b axis. All of the samples were enriched 
with 17 isotope by annealing in 17 02 gas at 900 C. We 
conducted most of the NMR measurements at 7 or 9 
Tesla. We assigned the three sets of 17 NMR signals 
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FIG. 1: Structure of two-leg ladder layer of A14CU24O41. 
Open and closed circles represent Cu and O atoms, respec- 
tively. 



to the 0(1) leg and 0(2) rung sites (see fig. |l|) and the 
additional oxygen site in the CUO2 chain layer based 
on the symmetry and temperature dependence of the 
v ' vq and NMR Knight shift tensors.el The quadrupole 
tensors, 17 ' 63 vq, are discussed below. 

In figure ||(a), we present the b-axis component of the 
oxygen nuclear quadrupole interaction, 17 t , Q[6], for the 
0(1) ladder leg site. Figure ||(b) shows the asymme- 
try parameter, 77, of Sri4Cu2404i as a function of tem- 
perature. A dramatic change in 17 vq occurs only for 
the doped samples, while the undoped LagCa8Cu2404i 
exhibits a minor decrease of 17 vq due to thermal ex- 
pansion. We define the temperature where this dra- 
matic change begins as T* and estimate T* = 210 K 
for Sri4Cu2404i and 140 K for SrnCa3Cu2404i. For 
Sr6CasCu2404i, 17 vq is changing even down to the low- 
est temperature (T* < 10 K). We also observed a similar 
dramatic change in temperature dependence at T* for 
17 vq at the 0(2) ladder rung site for these samples ijigure 
||) and for 63 vq at the copper site in Sri4Cu2404i.Qi3 Ex- 
tremely broad 63 Cu NMR lincshapcs in La6CagCu2404i 
and Sri4_ 2 ,Ca a; Cu2404i have prevented us from measur- 
ing 63 vq accurately. 

In order to understand the dramatic temperature de- 
pendence of 17 vq, we divide the electric field gradient, e<t. 
based on a standard ionic picture into two contributions,!?! 



eq = eq holc + (1 - 7)egi a 



(3) 



The first contribution, echoic, arises from holes in orbitals 
of the ion itself. An isotropic distribution of charge, such 
as from a filled electron shell, will not produce an electric 
field gradient, but unfilled shells can. If there is one hole 
in an oxygen 2p orbital, the electric field gradient at the 
nucleus is axial with the largest component along the lobe 
of the p orbital, vq = (^){^)e 2 Q < r~ 3 >. For our cal- 
culations, we take < r~ 3 >= 3.6|3.-ajtomic units, which i 



70% of the value for a free atom.BtL3 This results in a nu- 
clear quadrupole interaction of (2.66, -1.33, -1.33) MHz 
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FIG. 2: (a) Temperature dependence of 17 Vq(1)[6] at 0(1) 
ladder leg site in La6Ca8Cu2404i (°), Sri4Cu2404i (A), 
SrnCa 3 Cu2404i (□), and Sr 6 Ca 8 Cu2404i (•). (b) l7 r)(l) in 
Sri4Cu2404i. (Double r\ values caused by slight splitting in 
17 lines.) 
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FIG. 3: (a) Temperature dependence of 17 u Q (2)[b] at 0(2) 
ladder rung site in LaeCa8Cu2404i (°), Sri4Cu2404i (A), 
SniCa 3 Cu2404i (□), and Sr 6 Ca 8 Cu2404i (•). (b) 17 r?(2) in 

Sri4Cu24041- 
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for 17 2p x , for example. Similarly, one hole in the Cu 
30^2 _ y 2 orbital produces the nutdaar quadrupole interac- 
tion » (-38.5, -38.5, 77) MHz.Blia The second contribu- 
tion to the electric field gradient in equation (2), egiattice, 
arises from the charges of the other ions of the crystal. 
This electric field is then modified by the distortion which 
it creates in the electronic orbitals of the observed ion 
itself. This effect is accounted for by the Sternheimer 
antishielding factor, 7. For CuC>2 planes in high T c su- 
perconductors, the value of 7 is known to be 63 ^_Fp_. 20 
and 17 7 w -9 for 63 Cu and 17 0, respectivelyEM To 
calculate the lattice contribution to the electric field gra- 
dient, egiattice, we treat the ionic charges as point charges 
at their lattice locations.!! 3 ] The electric field gradient re- 
sulting from these point charges is then summed over the 
lattice out to a radius sufficient to achieve a stable value 
(typically 50 to 100 lattice spacings). For example, lattice 
point charge calculations of the electric field gradient for 
the undoped material, La6Ca§Cu2404i, were done with 
ionic charges for the (La,Ca) site +2.429, Cu +2, and O 
-2. Oxygen results for 7 = — 8 are (-193, +564, -371) 
kHz for the 0(1) site and (-713, +331, +381) kHz for the 
0(2) site, in rough agreement with the experimental val- 
ues (±60, ±516, +540±40) and (+690, ±170, ±540±40) 
kHz at 291 K. For the copper ladder site, using 7 = —20 
and an axial quadrupole inter.action from the one hole in 
the 3d a 2_ c 2 orbital of 77 MHzjLa point charge calculations 
give (+1.1, +16.9, -17.9) MHz in reasonable agreement 
with (+3.3, ±16.8, +13.5) MHz from measurements at 
room temperature. The oxygen quadrupole interaction, 
17 vq, of the undoped material, La6CasCu2404i, shows a 
steady decrease with increasing temperature. 

In contrast, for the doped materials, vq has a dra- 
matically different temperature dependence. In the fol- 
lowing paragraphs, we examine the possible causes for 
the temperature dependence and conclude that redistri- 
bution of the doped holes is responsible. In principle, 
the observed change in electric field gradient could arise 
from changes in one or more of the three parts of equa- 
tion (2), egi a tticej 7> and e^hoio- A change in the lattice 
contribution, egi a ttice 5 would be caused by a change in 
the crystal structure or lattice parameters as a function 
of temperature. As discussed by Carretta et al.,@ the 
change in lattice parameters as a function of temperature 
from 293 K to 520 K can account for less than 1% change 
in 63 vq, while the experimental change is about 15%. 
Similarly, we found that the change of lattice parameters 
would indicate a change of any component of 17 vq of 
1.5% or less, while the experimentally observed change 
is as much as 20%. In addition, the sign of the tem- 
perature dependence of 17 vq is opposite to what would 
be expected from lattice expansion. As the lattice ex- 
pands at higher temperatures, the electric field gradient 
from the lattice will generally decrease as we see in un- 
doped La6CasCu2404i. Clearly, the change in lattice 
parameters cannot account for the change in quadrupole 
interaction. Another option is that there is a local dis- 
tortion in the crystal structure. Our analysis using point 



charge calculations indicates that a local lattice distor- 
tion will not produce the experimentally measured sign 
of the temperature dependence at both the copper and 
two oxygen sites. The observed Cu G3 vq temperature de- 
pendence would require an increase in Cu-0 bond length, 
while the O 17 vq would require a decrease. More com- 
plicated lattice distortions involving further neighboring 
atoms do not overcome this. Thus, lattice changes are 
not responsible for the large change in electric field gra- 
dient. 

A second possible source for a change in electric field 
gradient is a change in thej-S|ternheimer antishielding fac- 
tor, 7. However, ShimizuEHl showed that 63 7 = —20 for 
copper does not vary significantly between copper oxide 
materials with doping and temperature. Furthermore, 
our point charge calculations reproduced both 63 vq and 
17 vq tensors of undoped La6CasCu2404i with 63 7 = —20 
and 17 7 = —8, very close to values known for high T c 
cuprates, 63 7 = —20, 17 7 = —9. Therefore, we conclude 
that temperature dependent 63 < 17 7 in hole-doped ladders 
is very unlikely. 

The third and only remaining possibility is that a 
change in the hole concentration in the ladder CU2O3 
layer changes eghoic of the doped samples. Unlike lattice 
changes, additional holes can account for the signs and 
symmetry of the 17 ' S3 vq temperature dependence. In the 
following, we will estimate the change in local hole con- 
centration in the CU2O3 ladder required to account for 
the changes in the electric field gradient. 

Additional holes affect the electric field gradient in two 
ways. First, a hole produces an electric field gradient at 
the nuclear site of its own atom through echoic of eq.(2). 
Second, a hole changes the charge of the ion, altering the 
charge environment of neighboring atoms. This changes 
the lattice electric field gradient, e<7i a tti C c7 for neighboring 
atoms. We calculate the effect of the holes on the neigh- 
boring atoms by point charge lattice summations using 
the values of 7 estimated for La6CasCu2404i. For our 
calculations, we consider the possibility of holes in the 
oxygen 2p orbitals of both oxygen sites and in the copper 
3d a 2_ c 2 or 3d 3b 2_ r 2 orbitals. The subscripts, a, b, and c, 
refer to the crystal axes where b is perpendicular to the 
ladder layer, a is along the rungs of the ladder, and c 
is along the ladder direction. We determine the amount 
of doped holes in each oxygen and copper orbital from 
63 vq and the tensor components of 17 vq(1,2) consider- 
ing consistently both the on-site effect of a hole on its 
own nucleus and the indirect effect of a hole on nearby 
atoms. The resulting equations relating the change in 
17 ' 63 vq to the change in hole concentration, Ah, in each 
orbital are (in units of MHz) 

A 17 ^ g (l)[6] = -1.33Afc 2po(1) - 1.33Aftap o(1 j 

- 0.55A/i Cu - 0.19Afr O (i) - 0.14A/i O (2) 

- 0.20A/i total (4a) 
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A 17 ^ Q (l)[c] = -1.33A/i 2pa(1) + 2.66A/i 2pc(1) 

+ 0.49A/i Cu + 0.12A/i o( i) + 0.05A/i O (2) 
+ O.lOA/itotai (4b) 

A 17 ^ Q (2)[fe] = -1.33A/i 2pa(2) - 1.33A/i 2pc(2) 

- 0.45A/i Cu - 0.29A/i O (i) - 0.Q4A/i O ( 2 ) 

- 0.17Afttotal (4c) 

A 17 ^ Q (2)[c] = -1.33A/i 2po(2) +2.66A/i 2po(2) 

- 0.30A/i Cu + 0.10A/i O (i) + 0.06Aft O ( 2 ) 
+ 0.10A/H o tai (4d) 



77A/i 3d 



2 _ c2 + 12Afe Cu + 35A/i (i) 
14Aft (2) + 13Afctotai (4e) 



A/lQ(l,2) — A/l 2pa ( lj2 ) + A/l2p c (l,2) 



(4f) 
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FIG. 4: Additional observed holes in the CU2O3 ladder 
layer per ladder Cu as deduced from 17 and 63 Cu vq for 
Sri4Cu2404i (a), SrnCa3Cu2404i (□), and Sr6CagCu2404i 
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A/icu = A/i 3da2 _ c2 (4g) 



A/itotai = (14/20)Aft C u + (14/20)Aft O (i) + (7/20)A/io(2) 

These equations have three types of terms. First, the 
terms that refer to specific atomic orbitals result from 
the on-site effect of a hole on its own nucleus. Second, 
the terms that refer to the holes on specific atomic sites 
result from the change in the lattice electric field gra- 
dient caused by holes on neighboring atoms in the lad- 
der layer. Third, the final A/i to tai term in these equa- 
tions results from the effect on the lattice electric field 
gradient of removing these holes from the chain layer. 
This third contribution is calculated based on the as- 
sumption that there are no holes in the ladder layer of 
Sri 4 Cu 2 404i in the low temperature region below 210 K 
and that additional holes in the ladder layer are removed 
evenly from all of the chain layer atoms. For the uniax- 
ially aligned powder samples, we have only measured b 
axis quadrupole frequency and we assumed that the dis- 
tribution of holes between the oxygen and copper orbitals 
was the same as in Sri4Cu 24 4 i. 

We summarize the total number of additional holes 
required in the ladder Cu 2 C>3 layer to produce the exper- 
imentally observed changes of 17 ' 63 vq in Figure [|, and 
the distribution of the holes between the oxygen 2p or- 
bitals and Cu 3d orbitals in table 1 for Sri4Cu 2 404i at 
500 K. The change in hole concentration for Sri4Cu 2 404i 
from low temperature to 500 K represents an observed in- 
crease of about one hole in the ladder layer out of the six 
holes implicit in the Sri4Cu 2 404i formula unit. Primar- 
ily, the holes go into the oxygen 2pcr orbitals, the 0(1) 
leg site 2p a (l) and 2p c (l) and the 0(2) rung site 2p Q (2). 



TABLE I: Additional observed holes in ladder layer orbitals 
of Sri 4 Cu240 4 i at 500 K deduced from 17 ' 6a v Q 



O 2p a (l) 


O 2p c (l) 


O 2p a (2) 


O 2p c (2) 


Cu 3d a 2_ c 2 


1.3±0.3 % 


1.3±0.3 % 


1.0±0.3 % 


0.3±0.3 % 


0.9±0.6 % 



There is possibly also some hole transfer to the ladder 
copper site, but there is much larger uncertainty in the 
calculation of the holes on the copper site because of the 
importance of the lattice contribution of the oxygen holes 
to the calculated 63 vq. 

We emphasize that figure ^ just shows the additional 
holes necessary to account for the temperature depen- 
dence of the electric field gradient. The low tempera- 
ture region of vq for Sri4Cu 2 404i was taken as a start- 
ing point from which the additional holes necessary to 
account for the temperature dependence of 17,63 ^q are 
calculated. This is done to avoid uncertainties in the lat- 
tice contribution to the electric field gradient. Lattice 
point charge calculations are clearly a simplified approx- 
imation which prevents an absolute determination of the 
hole concentration in the ladder Cu 2 03 layer. However, 
for SrgCa8Cu 2 404i, r 1 vq is still changing down to the 
lowest temperature (10 K) and our calculations do sug- 
gest that holes are still present in the ladder layer at 10 
K as shown in fig. [I| We emphasize that the physical 
picture given by this analysis is independent of the lat- 
tice point charge calculation and the value of 7. From 
Table 1, we see that the majority of the holes reside on 
the 0(1) site (2.6% on 0(1), 0.65% on 0(2), and 0.9% on 
Cu site per ladder Cu). For the 0(1) site, the on-site hole 
contribution is the dominant effect in the calculation and 
the contribution from the lattice effect of holes on nearby 
ions is < 25%. That the holes reside in oxygen 2pcr or- 
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bitals is consistent with x-ray absorption spectroscopy^ 
and other copper, o^ide materials, such as YBa2Cu307 
and (La,Sr)Cu0 4 Hi2 



III. DISCUSSION 

Other measurements of these A14CU24O41 materials 
also point towards changes in hole concentration at T*. 
Optical conductivity measurements^ of the Drude peak 
at low frequency are attributed to carriers in the lad- 
der layer. The temperature dependence of the integral of 
that peak, which represents N/m* (carrier number over 
effective mass) , is similai^a.that of the electric field gra- 
dient. Charge transportE3EZl show an anomaly at around 
T* for several doping levels as summarized in figure ^ 
and also a-possible collective excitation in Sr^Gug^X^i 
below T*x3 In addition, magnetic susceptibilityllSO in- 
dicates that the ladder layer contributes to the suscep- 
tibility only above T*. The mean free path of ladder 
magnetic excitations also shows a sharp decrease around 

jl9| 

Also correlated with the changes in electric field gradi- 
ent are magnetic changes seen in .the nuclear spin lattice 
relaxation rate, 1/T l5 for 63 Cu.B 1/Ti generally mea- 
sures low energy spin fluctuations. As shown in figure 
^(b), 63 Cu 1/Ti crosses over from a low temperature 
gapped regime where 1/Ti roughly follows an activa- 
tion law, exp(— Ati/ZcbT), to a paramagnetic regime. 
The onset of the change in the electric field gradient, 
T*, is at the beginning of the crossover from the low 
temperature gapped regime to the paramagnetic regime. 
In addition, for Sri4Cu240 4 i, Takigawa et alJ3 and Car- 



FIG. 6: (a) Temperature dependence of the ratio 
M (l/Ti)/ a3 (l/Ti) for Sr 14 Cu 24 04i (A), SmCasCu^CUi (□), 
and Sr6Ca8Cu24C>4i (•). A few representative error bars 
shown, (b) 63 (1/Ti) for La 6 Ca 8 Cu240 4 i (o), Sri 4 Cu 2 40 4 i 
(A), SrnCa3Cu2404i (□), and Sr6CagCu2404i (•). Lines are 
guides for the eye. 

retta et all noted that for T < 225 K, the fit of the 
nuclear magnetization decay to the standard solution of 
the rate equations becomes poor, and the ratio ( 65 1/Ti) 
I decreases from ( e5 j/ 63 j) 2 = 1.15 indicating 

that charge fluctuations are contributing to the relax- 
ation. The quadrupolar relaxation is maximum at about 
100 K for Sri4Cu2404i indicating that charge fluctua- 
tions have slowed down to the NMR frequency.EH3 As 
shown in Figure |^(a) , we found that quadrupolar relax- 
ation is significant below T An for all hole doped 
samples. Most likely, these charge fluctuations are asso- 
ciated with very slow hole motion in the ladder plane. 

The additional holes that are observed in the CU2O3 
ladder above T* could arise from two obvious sources: 
transfer from the chain layer, or derealization from un- 
observed regions of the ladder layer. Transfer of holes 
from the chain layer seems quite plausible because of 
the coincidence of the temperature of the derealization 
of holes on the chain with the appgarance-of the holes 
in the ladder layer at J 1 *. NMR,Ll ESR,EJ and x-ray 
diffraction experimentsEil all show that the holes in the 
chain delocalize at around T* = 210 K for Sri4Cu2404i. 
The ESR experiment also indicates charge ordering of 
the chain at 170 K fpc Sri2Ca2Cu2404i and at ~ 80 
K for SrgCa5Cu2404i,E3 consistent with our observation 
of hole increase in the ladder layer starting at 140 K 
for SrnCa3Cu2404i. Transfer of holes from the chain 
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to the ladder layer also explains one apparent discrep- 
ancy. n Below the chain chacge ordering temperature, 
NMR0 and neutron scatteringc3 measurements show sep- 
arated dimers on the chain of Sri4Cu2404i. The sepa- 
rated dimer model has 6 holes per formula unit, which is 
all of the holes expected in Sri4Cu2404i. If no holes 
move between chain and ladder layers, this does not 
leave any holes for the ladder layer. This conflicts with 
the measurement of roughly 1 hole per formula unit in 
the ladder plane of Sr 14 Cu2404i at high temperature tw 
this work and optical conductivityu Takigawa, et alu 
noted this discrepancy and also pointed out that any non- 
stoichiometry did not seem large enough to account for 
this effect. On the other hand, if holes do move from the 
chain layer to the ladder layer with increasing tempera- 
ture, there is a consistent picture: the holes are all on the 
chain below 210 K and gradually transfer to the ladder 
plane at higher temperatures, as shown in Fig. ||. 

The second possible source of holes is derealization 
from unobserved regions of the ladder layer. The hole 
concentration in the CU2O3 ladder below T* can be 
smaller for the majority of the Cu and 0(1,2) sites, if 
holes are localized within the NMR time scale (~ /xsec) 
below T*. One possible piece of evidence in support 
of this possibility is the existence of a collective exci- 
tation in the conductivity of Sri4Cu2404i below ^170 
K.L3 This might indicate the motion of pinned charges 
in the ladder layer. However, the origin of the collec- 
tive excitation is not clear. In addition, the spectral 
weight of the conductivity peak is at least six orders of 
magnitude smaller than that expected using Pl ~ 0.07 
and an effective mass equal to the free electron mass. 
Thus, this excitation may involve a very small number 
of holes. The anomalies in charge transport and mag- 
netic susceptibility measurements have been interpreted 
as charge localization at T*J13 but this does not seem 
to account for the NMR results. This scenario requires 
NMR line broadening or splitting around T* due to the 



localization of holes, but we did not observe any signifi- 
cant broadening of the NMR lineshape around T*. Some 
NMR quadrupole satellite lines are slightly split, but the 
size of the splitting is temperature independent (shown 
by X) in figure ^) . Additionally, the peak in quadrupole 
relaxation for Sri4Cu2404iLj implies that the hole motion 
slows down to the NMR frequency only at around 100 K, 
so the charge dynamics at higher temperatures such as 
T*=210 K should be faster than the NMR time scale, 
not slower. Therefore, this derealization scenario does 
not seem as consistent as the possibility of holes trans- 
ferring from the chain to the ladder layer. However, the 
quadrupolar relaxation does imply that the ladder layer 
has nearby charge fluctuations, so it is possible that some 
holes do localize on the ladder layer. 

IV. CONCLUSION 

We report a detailed study of 17 and 63 Cu vq mea- 
surements on A14CU24O41 two-leg ladder sites. The 
charge environment of the ladder layer for the doped sam- 
ples dramatically changes at temperatures above T*. An 
increase in the effective hole concentration is observed in 
the ladder layer above T*, primarily in oxygen 2pa or- 
bitals. We suggest that these holes may be transferred 
from the chain layer. The change in the charge environ- 
ment of the ladder layer occurs at the beginning of the 
magnetic crossover from the low temperature gapped re- 
gion to the high temperature paramagnetic region. This 
suggests that spin and charge degrees of freedom are 
closely connected in the hole doped two-leg CU2O3 ladder 
of A14CU24O41. Holes are excluded from the majority of 
the ladder plane below T*, when the two-leg ladder is in 
the low temperature spin-gapped state. Above T*, the 
spin-gapped state is disturbed by holes and spin excita- 
tions are seen in 1 /Ti and magnetic susceptibility. 

Early parts of this work were supported by NSF DMR 
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